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_ , ' Abstract 

> 

^^ • The prospects of measuring the leptonic angles and CP-odd phases at a neutrino 

factory are discussed in the scenario of three active plus one sterile neutrino. We 

^^ I consider the v^ — > u^ LSND signal. Its associated large mass difference leads to 

observable neutrino oscillations at short (-^ 1 km) baseline experiments. Sensitivities 

Qv . to the leptonic angles down to 10~^ can be easily achieved with a 1 Ton detector. 

"^ I Longer baseline experiments (~ 100 km) with a 1 Kton detector can provide very 



p i' clean tests of CP-violation especially through tau lepton detection. 
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d • 1 Introduction 



Indications in favour of neutrino oscillations have been obtained in both the at- 
mospheric and solar neutrino experiments [1,2,3,4]. The atmospheric neutrino data 
require Ar/i^^^ ~ (0.3 — 7) x 10~^ eV^ whereas the solar neutrino data prefer either 
Aml^^ ~ 10~^ — 10~^ eV^ if interpreted as MSW (matter enhanced) oscillations, 
or ~ 10"^° eV^ if the vacuum solution is preferred by nature. The LSND data 
[5] would indicate a z/^ — *> z/g oscillation with a third, very distinct, neutrino mass 
difference: Am^^^yD ~ 0-3 — 1 eV^ . The usual three neutrino picture has then to 
be enlarged to explain the present ensemble of data: four different light neutrino 
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species are needed. The new light neutrino is usually denoted as sterile [6,7], in it 
cannot have the standard weak interaction charge in order to comply with the strong 
bounds on the Z° invisible decay width [8]. 

We derive here convenient parametrizations of the physical mixing angles and CP 
phases in the case of four neutrino species, and study their experimental signals at 
the neutrino factor-ij^^ 

For the sake of illustration, we shall consider as a "reference set-up" the neutrino 
beams resulting from the decay of n^ = 2 x 10^°/i"'"s and/or fi~s in a straight 
section of an E^ = 10 — 50 GeV muon accumulator ring. Muon energies of 40 — 50 
GeV are at present under discussion [14] as a convenient goal, as they allow good 
background rejection [15] and do not preclude the exploration of neutrino signals at 
lower energies. This is because the number of neutrinos in a given energy bin does 
not depend on the parent muon energy, while the total number of oscillated and 
interacted neutrinos increases with E^ [16]. As the dominant signals are expected 
to peak at "LSND" distances, most of the parameter space can be explored in 
experiments at short baseline (SBL) distances of 0(1) km and with a very small 
1 Ton detector. To optimize the observability of CP-violation asymmetries larger 
distances (0(10-100) km) and larger detector (1 kTon) are needed. 

In section 2 we describe our parametrization. Sections 3 and 4 discuss the experi- 
mental sensitivities to angles and CP phases, respectively. We conclude in section 
5. 



2 Four-Neutrino species 



The masses of the four neutrino species can be ordered in two ways, as depicted in 
Fig.l: 

(1) Three almost degenerate neutrinos, accounting for solar and atmospheric oscil- 
lations, separated from the fourth neutrino specie by the LSND mass difference; 

(2) Two almost degenerate neutrino pairs, accounting each for solar and atmo- 
spheric oscillations, separated by the LSND mass gap. 

We refer to these possibilities respectively as class-I and class-II schemes. Both of 
them contain a number of sub-classes accounting for all the possible ordering of 
the neutrino species. It has been shown in [17] that the combined analysis of solar, 
atmospheric and LSND experiments tends to exclude all class-I schemes. We consider 
then two possible class-II schemes and choose for definiteness the hierarchycal class- 
IIB scenario in Fig.l. The lighter pair (1-2) is separated by the solar mass difference 



For a detailed analysis of the three family scenario see [9,10,11,12,13]. 
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Fig. 1. Different type of four-neutrino families scenarios: Class-I scenarios (left); Class-II 
scenarios (right). 

and the heavier pair (3-4) by the atmospheric oneQ. We work in the convention: 

I/, ~ i/i ^ z/e ~ V2 Am^2 = ^"^L (1) 



with the solar deficit assigned to v^ — z/g oscillations, and 

Vt^^V^^Vr^ va Am^4 = Am^tm 



(2) 



with the atmospheric anomaly due to z^^ — Vr oscillations. The alternative possibility 
of identifying the atmospheric anomaly as a z/^ — Vg oscillation is actually disfavoured 
by SK data[l]. 

Whatever the mechanism responsible for neutrino masses, given n light neutrino 
species oscillation experiments are only sensitive to a unitary n x n mixing matrix 
"a la CKM" . In all generality, the parameter space of a four-species scenario would 
consist of six rotation angles and three complex phases if the neutrinos are Dirac 
fermions, while it spans six angles plus six phases if the neutrinos are Majorana 
fermions. Among the latter, three are pure Majorana phases and thus to be disre- 
garded in what follows as they are unobservable in oscillations, reducing the analysis 
to the mentioned 4x4 "Dirac-type" system. 

Within our conventions and the large mass hierarchy indicated by data: 



AmL 



Am?2 < Aml^ = Amj^ < Am 



2 
LSND 



Am23, 



(3) 



^ The alternative option, with a small separation in the heavier pair and a larger one for 
the lighter pair, amounts to reversing the sign of the LSND mass difference. The analysis 
below does not depend on this choice. 
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Table 1 

Parameter space for four neutrino families: for Dirac neutrinos we consider the general 
case with three non-zero mass differences and the particular case (considered in the rest 
of the paper) with one or two mass differences set to zero; for Majorana neutrinos we 
consider only the general case. 

two approximations are useful: 



(1) Am^2 = Am34 = 0, "one mass scale dominance" (or minimal) scheme, 

(2) Am^2 = 0) "two mass scale dominance" (or next-to-minimal) scheme. 



The number of independent angles and phases is then reduced as reported in Tab. 1. 
The minimal scheme is sufficient to illustrate the sensitivity to the mixing angles 
other than the solar and atmospheric ones. The next-to-minimal scheme is necessary 
to address the question of CP-violation, as two non-zero mass differences are required 
to produce observable effects, alike to the standard three-family scenario. 

A general rotation in a four dimensional space can be obtained by performing six 
different rotations f/^ in the (i,j) plane, resulting in plenty of different possible 
parametrizations of the mixing matrix, disregarding phases. We choose the following 
convenient parametrization, given the hierarchy of mass differences in eq. (3): 



U = f/l4(^14)f/l3(^13)f/24(^24)f/23(^23,53)f/34(^34'52)f/l2(^12,'5l) 



(4) 



As shown in Table 2, if a given mass difference vanishes the number of physical 
angles and phases gets reduced by one. A convenient parametrization of the angles 
is that in which the rotation matrices corresponding to the most degenerate pairs of 
eigenstates are placed to the extreme right. If the eigenstates i and j are degenerate 
and the matrix Uij is located to the right in eq. (4), the angle dij becomes automat- 
ically unphysical. When a different ordering is taken no angle disappears from the 
oscillation probabilities. A redefinition of the rest of the parameters would then be 
necessary in order to illustrate the remaining reduced parameter space in a trans- 
parent way. Our parametrization corresponds thus to the "cleanest" choice, having 
settled at the extreme right the rotations corresponding to the most degenerate 
pairs. 



In the "one mass dominance" scheme, the pairs {Ou, ^i) and (6*34, ^2) decouple auto- 
matically. In the "two mass dominance" scheme, only the pair (^12,(^1) does. Thus 
only the exact number of physical parameters, according to Table 2, remains both 
in the minimal and next-to-minimal schemes. Notice that it is also important to 
distribute the phases so that they decouple, together with the angles, when they 
should. 



2. 1 Sensitivity reach of the neutrino factory for four neutrino species 



We concentrate now on the sensitivity to the different angles in the "one mass scale" 
approximation, discussed at the beginning of this section. Four rotation angles (6'i3, 
^14, ^23 and 624) and one complex phase {5^) remain. The two rotation angles that 
have become unphysical are already tested at solar [612 in our parametrization) 
and atmospheric (^34) neutrino experiments. The remaining four can be studied 
at the neutrino factory with high precision, due to the rich flavour content of the 
neutrino beam. Notice that the number of measurable flavour transitions is enough 
to constraint them. Consider the following channels: 



'^e ^ t'^ ^ Ai^ (yU^appearance) 
i/^ ^ z/^ ^ ji^ (/i~ disappearance) 



V^^ T 



(reappearance) 



u^^Ut-^t (r appearance). (5) 

Separate the CP-even terms, P^p , from the CP-odd ones, Pp^ , in the following way: 

P{Ua -^ Vp) = Pcp{l^a -* l^fi) + Pofi^^a ^ l^p)- (6) 

The CP-even parts acquire simple forms: 



Pcp{iye ^ jy^) =^0^3^24^23^23 S"^ ( ^^ I > iV 

Pcp{iy^ -* z/^) = 1 - 444(s^3 + 5^3^) sin^ -^^ J , (8) 
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(9) 



(10) 



Notice that the physical phase appears in Pcpii^e -^ t-'t) and Pcp{i^tJ. -^ ^r) in a 
pure cosine dependence. Actually, no CP-odd observable can be build out of the 
oscillation probabilities in this approximation in spite of the existence of a physical 
CP-odd phase in the mixing matrix. 

The existing experimental data impose some constraints on the parameter space, 
although the allowed region is still quite large. Bugey and Chooz are sensitive to the 
oscillation z/p ^ z/p. 






P(Ue ^Ue) = l- ^cl.c\As\. + S^, C^ ) siu^ ^-^ 1 , (ll^ 



V 4E 



resulting in the bound 



c^3 sin^ 2^24 + C24 sin^ 2^23) < 0.2 , (12) 



while Bugey gives a slightly stronger constraint in the larger mass range allowed by 
LSND. In our computations we safely stay within both experimental constraints. 
Notice that in the assumption of small angles, this bound forces the mixings S23 and 

S24 to be small and leaves more freedom in the mixings of the sterile neutrino: s^3 

2 

14- 



and s^ 



The allowed mass range for the LSND signal of z/g -^ z/^ transitions provides the 
constraint 10~^ < c\^c^^ sin^ 2^23 < 10~^, which fits nicely with the Chooz constraint 
to point towards small S23 values. 

We choose to be "conservative", or even "pessimistic", in order to illustrate the 
potential of the neutrino factory. In the numerical computations below we will make 
the assumption that all angles crossing the large LSND gap, 6'i3, ^14, ^23 and 6*24 are 
small. 

The large LSND mass difference, Am23 — leV^, calls for a SBL experiment rather 
than a long baseline one. Consider for illustration a 1 Ton detector located at ~ 1 
km distance from the neutrino source. We consider a muon beam of Ey^ = 20 GeV, 
resulting in Ncc — 10^ charged leptons detected, for a beam intensity of 2 x 10^° 
useful fi~ per year. An efficiency of e^ = 0.5, e,- = 0.35 for /i and r detection 
respectively, and a background contamination at the level of 10~^Ncc events are 
included. 

• ^23 and ^13 from fi channels 

The n~^ appearance channel is particularly sensitive to ^23- Fig. 2 shows the sensitiv- 
ity reach in the S23/Am23 plane for different values of ^13. Inside the LSND allowed 
region the dependence on ^13 is mild: S23 can reach 10~^ for ^13 ~ 1° or 6 x 10^® for 
^13 - 60°. 




Fig. 2. Sensitivity reach in the s^s/ATTigs plane at different values of 6i^ = 1°,10°,30° 
and 60° for fj.'^ appearance. We consider a 1 Ton detector at 1 km from the source and 
2 X lO^*' useful muons/year. 

Concerning the sensitivity to 613, Fig. 3 (left) illustrates the reach from fi appearance 
measurements in the sl^/Arn^^ plane, for different values of 623- Inside the LSND 
allowed region, the sensitivity to this angle strongly depends on the value of 6^23, 
with the larger sensitivity attained for large values of 623, a scenario somewhat 
disfavoured by the LSND measurement. For small values of 623 (— 1°), the smallest 
testable value of 5^3 is ~ 10~^. Nevertheless, in this range the muon disappearance 
channel proves quite more sensitive: in Fig. 3 (right) the sensitivity goes down to 
values of 5^3 as small as 10~^ for 623 — 1°. 

• 614 and 62A from r channels 

The T~ appearance channel is quite sensitive to both 5^4 and S24. Fig. 4 illustrates 
the sensitivity to 5^4 as a function of 6*13: for about 1°, sensitivities of the order 
of 10~^ are attainable, while for 10° the reach extends to 4 x 10~^. For even larger 
values of 6'i3 it goes down to 10~^ (we recall that 6*13 is not constrained by the present 
experimental bounds). 



Fig. 5 (left) depicts the foreseeable sensitivity reach to S24 as a function of 6^23 
small values of 6^23 the sensitivity to S24 goes down to 10~^. 



for 



In contrast, the r"*" appearance channel looks less promising. This is illustrated in 
Fig. 5 (right): due to the relative negative sign between the two terms in the analytic 
expression for P{i'e — *■ i^t), eq. (9), cancellations for particular values of the angles 
occur, resulting in a decreasing sensitivity in specific regions of the parameter space. 
For instance, for 6'23 = 10°, the reach in 534 splits into two separate regions for the 
LSND allowed range Ar?223 ~ 10^ eV^. This sensitivity suppression is absent in the 
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Fig. 3. Sensitivity reach in the Si^/Am2^ plane at different values 0/623 = 1°,10° and 
30° for n^ appearance (left) and disappearance (right). We consider a 1 Ton detector at 
1 km from the source and 2 x 10^'^ useful muons/year. 

10' 



< 



',,= 30" 



)„= lO" 




E„= 20 GeV 



Size = 1 Ton 



10"' lO"" 10"' lO"- 10"' 10"' 10"' 10" 

sin^e,. 



Fig. 4. Sensitivity reach in the Si^/Am23 plane at different values of 9 13 = 1°,10° and 
30° for T^ appearance. We consider a 1 Ton detector at 1 km from the source and 2 x 10^'' 
useful muons/year. 

T^ channel as the relative sign between the two terms in P(z/^ -^ Vt)-, eq. (10), is 
positive. 

The overall conclusion of this analysis is that, in the minimal scheme for four- 
neutrino families, a 1 Ton near detector with /i and r charge identification is suitable 
to fully explore the CP-even part of the parameter space. 
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Fig. 5. Sensitivity reach in the s^Jt^m^.^ ■plane at different values 0/623 = 1°,10° and 
30° for T~ (left) and r"*" (right) appearance. We consider a 1 Ton detector at 1 km from 
the source and 2 x 10^'^ useful muons/year. 

2.2 CP Violation with four light neutrino species 



As in the standard three-family scenario, in order to reach observable CP-odd effects 
in oscillations it is necessary to have both physical CP-odd phases and at least two 
non- vanishing mass differences. The next-to minimal or "two mass scale dominance" 
scheme, described at the beginning of this section, is thus suitable. 

As explained above, the parameter space consists of five angles and two CP-odd 
phases. Expanding the transition probabilities to leading order in Am^^^ (i.e. Am^^ 
in our parametrization) , it follows that their CP-odd components are^ : 
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^ At this order also sub-leadings in the CP-even sector contribute. Although we do no 
illustrate them, all orders in Atti^^^^ are included in the numerical computations. 






Two distinct phases appear, 62 and 63, in a typical sinusoidal dependence which is 
the trademark of CP-violation and ensures different transition rates for neutrinos 
and anti neutrinos. 

CP-odd effects are observable in "appearance" channels, while "disappearance" ones 
are only sensitive to the CP-even part. The latter is mandated by CPT [18]. In 
contrast with the three- neutrino case, the solar suppression (see [9]) is now replaced 
by a less severe atmospheric suppression. CP-violating effects are then expected to 
be one or two orders of magnitude larger that in the standard case, and independent 
of the solar parameters. 

Staying in the "conservative" assumption of small ^13,^14,^23)^24) we compare two 
democratic scenarios, in which all these angles are taken to be small and of the same 
order: 

(1) Set 1: ^34 = 45°, % = 5° and Am^^^ = 2.8 x IQ-^ eV^ for Am^^jv^ = 0.3 eV^; 

(2) Set 2: ^34 = 45°, % = 2° and Am^^^ = 2.8 x IQ-^ eV^ for Am^^jv^ = 1 eV^. 

The value chosen for Am^^^ is the central one of the most recent SuperK analysis 
[1]. For illustration we consider in what follows a 1 kTon detector located at 0(10) 
km distance from the neutrino source. In the figures below, the exact formulae for 
the probabilities have been used. 

The easiest way to measure CP-violation in oscillation is to built a CP-asymmetry 
or a T-asymmetry: 
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-^ ^a) 
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A^p and A^^ are theoretically equivalent in vacuum due to CPT, and matter effects 
are negligible at the short distances under consideration. Their extraction from data 
at a neutrino factory is quite different, though. Consider, as an example, the (e/i)- 
channel. The CP-asymmetry, A^, would be measured by first extracting P{i'^ -^ 
z/g) from the produced (wrong-sign) /i~s in a beam from /i^ decay and P{ue -^ ^n) 
from the charge conjugate beam and process. Notice that even if the fiuxes are 
very well known, this requires a good knowledge of the cross section ratio (t(z/^ — > 
ix^)/a{u^ — > /i~). Conversely, the measurement of the T-asymmetry, Aj^, requires 
to consider P(z/^ — > u^) and thus a good e charge identification, that seems harder 
to achieve. In the following we will deal only with CP-asymmetries. 
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Fig. 6. Signal over statistical uncertainty for CP violation, in the v^ — > v^ channel, for 
the two sets of parameters described in the text (Set 1 on the left and Set 2 on the right). 
We consider a 1 kTon detector and 2 x lO^*' useful muons/year. 

A central question on the observability of CP-violation is that of statistics. We do 
not exploit here the explicit E^ dependence of the CP-odd effect, and we consider 
the neutrino-energy integrated quantity: 



<"(5) 



{iV[/i-]/iV,[e-]}+-{iV[/i+]/iVje+]}_ 
{iV[;.-]/iV4e-]}+ + {iV[;.+]/iV4e+]}_ 



(19) 



where the sign of the decaying muons is indicated by a subindex, N[^'^\ {N[^~\) are 
the measured number of wrong-sign muons, and A''o[e"'"] (A^o[e~]) are the expected 
number of ^^{ve) charged current interactions in the absence of oscillations. In order 
to quantify the significance of the signal, we compare the value of the integrated 
asymmetry with its error, in which we include the statistical error and a conservative 
background estimate at the level of 10~^. 

The size of the CP-asymmetries is very different for /i channels and r channels. For 
instance for Set 2, they turn out to be small in z/g — v^j, oscillations, ranging from the 
per mil level to a few percent. In contrast, in z/^ — Vr oscillations they attain much 
larger values of about 50% — 90%. This means that their hypothetical measurement 
should be rather insensitive to systematic effects, and other conventional neutrino 
beams from pion and kaon decay could be appropriate for their study. 

• /i appearance channels 

Fig. 6 shows the signal over noise ratio for the integrated CP asymmetry, eqs. (19), 
in the wrong sign muon channel, that is i^e -^ ^fi versus Ue —>■ v^ oscillations, as a 
function of the distance. Matter effects, although negligible, have been included. For 
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Fig. 7. Signal over statistical uncertainty for CP violation, in the v^ — > Vr channel, for 
the two sets of parameters described in the text (Set 1 on the left and Set 2 on the right). 
We consider a 1 kTon detector and 2 x 10^^ useful niuons/year. 

the scenario and distances discussed here, the scahng laws are analogous to those 
derived for three neutrino species in vacuum, [9,10,11], that is 



aCP 

^-R— oc ^|E 



e/i 



sm 



4K. 



(20) 



The maxima of the curves move towards larger distances when the energy of the 
muon beam is increased, or the assumed LSND mass difference is decreased. More- 
over, increasing the energy enhances the significance of the effect at the maximum 
as expected. At E^ = 50 GeV, 6 standard deviation (sd) signals are attainable at 
around 100 km for the values in Set 1, and just 2.5 sd at 30 km for Set 2, levelling off 
at larger distances and finally diminishing when E^/L approaches the atmospheric 
range. 



u-r versus u, 



M 



oscillations 



• r appearance channels 

In Fig. 7 we show the signal over noise ratio in z/^ 
as a function of the distance. The experimental asymmetry and the corresponding 
scaling law are obtained from eq. (19) and eq. (20), with the obvious replacements 
e ^ /i and /i -^ r. A larger enhancement takes place in this channel as compared 
to the Ue — > i'^ one: over 60 sd for Set 1 and 33 sd for Set 2 are a priori attainable. 
These larger factors follow from the fact that the CP-even transition probability 
Pcpi^fi^r) is smaller than Pcp{i^e^fi), due to a stronger suppression in small mixing 
angles. Notice that the opposite happens in the 3-species case. Bilenky et al. [17] 
had previously concluded that the r channel was best for CP-studies in the four 
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Fig. 8. CP violation asymmetry in the I'e -^ i^^i Oeft) and v^ -^ v^- (right) channel for 
E^ = 20 Ge V, angles and mass differences as in Set 2 and for different choice of the CP 
phases: Si = 62 = S^ = 7r/2 (full line), 7r/4 (dashed line) and tt/12 (dotted line). We 
consider a 1 kTon detector from the source of a 2 x 10^^ muon/year beam. 

species scenario. Their argument relied, though, on the fact that the parameter 
space involved in Ur oscillations is experimentally less constrained than the z/^ one, 
a freedom we have not used here, staying within the more natural assumption that 
all the angles in the next-to-minimal scheme, except the atmospheric one, 634, are 
small. 



The results in the Uf. — > u^- channels are almost identical to the z/g 
deserving a separate discussion. 



Un ones, not 



The phase dependence is shown in Fig. 8, with the expected depletion of the signal 
for small CP phases. For small values of the phases, i.e. 61 = 62 = S3 = 15°, the 
significance drops to the la level. 



3 Conclusions 



The ensemble of solar, atmospheric and LSND neutrino data, are analysed in a three 
active plus one sterile scenario. A neutrino factory from muon storage beams has 
much higher precision and discovery potential than any other planned facility. The 
reach of SBL experiments is extremely large. We have derived one and two mass 
scale dominance approximations, appropriate for CP-even and CP-odd observables, 
respectively. The number of useful observables is sufficient to determine or very 
significantly constrain all the mixing parameters of a four-generation mixing scheme. 



13 



CP violation may be easily at reach, specially through "r appearance" signals. In 
these channels the CP-asymmetries are so large that even neutrino beams from 
conventional pion and kaon decays may be sufficient for their detection. Matter 
effects are irrelevant in the energy integrated CP-odd observables for such a short 
baseline, making the measure of CP-violation extremely clean in this four neutrino 
scenario. 
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